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Abstract 

We construct a family of second-order linear difference equations parametrized by the hypergeomet- 
ric solution of the elliptic Painleve equation (or higher-order analogues), and admitting a large family of 
monodromy-preserving deformations. The solutions are certain semiclassical biorthogonal functions (and 
their Cauchy transforms), biorthogonal with respect to higher-order analogues of Spiridonov's elliptic beta 
integral. 
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H ■ 1 Introduction 



In |20j , Sakai introduced an elliptic analogue of the Painleve equations, including all of the known discrete 
(and continuous) Painleve equations as special cases. Unfortunately, although Sakai's construction is quite 
natural and geometric, it does not reflect the most important role of the ordinary Painleve equations, namely 
as parameters controlling monodromy-preserving deformations. 

As with the ordinary Painleve equations, the elliptic Painleve equation admits a special class of "hyperge- 
ometric" solutions [HI [17] that in the most general case can be expressed via n-dimensional contour integrals 
with integrands expressed in terms of elliptic gamma functions. It is thus natural, as a first step in constructing 
an isomonodromy interpretation of the elliptic Painleve equation, to attempt to understand that interpretation 
in the hypergeometric case. 

In the present work, we do precisely that: associated to each elliptic hypergeometric solution of the elliptic 
Painleve equation, we construct a corresponding second-order linear difference equation that admits a family of 



discrete "monodromy-preserving" deformations. (In fact, the construction works equally well for higher-order 
analogues of the relevant elliptic hypergeometric integrals.) The construction is based on an analogue of the 
approach in [HIITT]. There, a linear differential equation deformed by the hypergeometric case of the Painleve 
VI equation is constructed as a differential equation satisfied by a family of "semiclassical" (bi-)orthogonal 
polynomials. Our construction is much the same, although there are several technical issues to overcome. 

The first such issue is, simply put, to understand precisely what it means for a deformation of an elliptic 
difference equation to preserve monodromy, or even what the monodromy of an elliptic difference equation is. 
While we give only a partial answer to this question, we do define (in Section 2 below; note that many of 
the considerations there turn out to have been anticipated by Etingof in [8]) a weakened form of monodromy 
that, while somewhat weaker than the analogous notions at the g-difference [9] and lower [4] levels, is still 
strong enough to give a reasonably rigid notion of isomonodromy deformation. Indeed, two elliptic difference 
equations have the same weak monodromy iff the corresponding difference modules (see |16| ) are isomorphic; 
the same holds for ordinary difference equations, even relative to the stronger notion of monodromy [5] . The 
key observation is that a fundamental matrix for a p-elliptic g-difference equation is also a fundamental matrix 
for a g-elliptic p-difference equation; this latter equation (up to a certain equivalence relation) plays the role of 
the monodromy. (The result is similar to the notion of monodromy introduced by Krichever in [13j ; while our 
notion is weakened by an equivalence relation, it avoids any assumptions of genericity.) 

In Section 3, we develop the theory of semiclassical elliptic biorthogonal functions, functions biorthogonal 
with respect to a density generalizing Spiridonov's elliptic beta integral [22] by adding m additional pairs of 
parameters. The key observation is that such functions can be constructed as higher-order elliptic Selberg 
integrals of a special form; in addition, their "Cauchy transforms" can also be so written. This gives rise to 
several nice relations between these functions, which we describe. Most important for our purposes is their 
behavior under p-shifts; the biorthogonal functions themselves are p-elliptic, but if we include the Cauchy 
transforms, the overall action is triangular. We can thus construct from these functions a 2 x 2 matrix which 
satisfies a triangular g-elliptic p-difference equation, analogous to the Riemann-Hilbert problem associated to 
orthogonal polynomials ( |10|, §3.4]; see also [6] for a general exposition). By the theory of Section 2, this 
immediately gives rise to a p-elliptic g-difference equation, and symmetries of the p-difference equation induce 
monodromy-preserving deformations of the g-difference equation. 

Finally, in Section 4, we compute this difference equation and the associated deformations. Although we 
cannot give a closed form expression for the difference equation, we are able at least to determine precisely 
where the difference equation is singular, and at each such point, compute the value (or residue, as appropriate) 
of the shift matrix. Together with the fact that the coefficients are meromorphic p-theta functions, this data 
suffices to (over)determine the shift matrix. 

In a foUowup paper [3] , with Arinkin and Borodin, we will complete the isomondromy interpretation of the 
elliptic Painleve equation by applying the ideas of [2] to show that any difference equation having the same 
structure as the ones constructed below admits a corresponding family of monodromy-preserving deformations, 
and moreover that (when m — I) Sakai's rational surface can be recovered as a moduli space of such difference 
equations. 
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Notation 

The elliptic Gamma function [TH] is defined for complex numbers p, q, z with |p|, \q\ < 1, z 7^ 0, by 

^vA^^ ■■- 11 1 _„.„,. ^ (1-1) 



1 — p'q^ z 



and satisfies the relations 



where 



satisfies 



rp,,(pz) = B^{z%^^{z) (1.2) 

I^,g((7z) = 6p{z)Vp^q{z) (1.3) 

r,,,(pg/z) = rp,,(z)-i (1.4) 

IJ),9(M2)IJ,g(^) = -z"^IJ,,q(pz)i;,,(gz), (1.5) 

^pW:-n(l-P'^V^)(l-p'^) (1-6) 

0<i 

0p(z) = -z0p(l/z) (1.7) 

- 0p(p/z). (1.8) 

By convention, multiple arguments to a Gamma or theta function represent a product; thus, for instance 

IJ,,(moz='=^) = I^,9(uoz)I^,g(uo/z). (1.9) 

We will also make brief use of the third-order elliptic Gamma function 

^UM--= n (i-pV*'^)(i-p'^V^'*'^V^), (1-10) 

0<ij\A: 

which satisfies 

r+,,*(te) = rp,,(x)r+^,,(x), (i.ii) 

Yl^^,(pqtlx)=Vl^^,{x)- (1.12) 

for our purposes, this appears only as a normalization factor relating the order 1 elliptic Selberg integral to the 
hypergeometric tau function for elliptic Painleve. 



2 Elliptic difference equations 

Let p be a complex number with |p| < 1. A (merom.orphic) p-theta function of multiplier az'^ is a meromorphic 
function f{z) on C* := C \ {0} with the periodicity property ,f{pz) = az^ f{z). (To justify this definition, 
observe that the composition /(exp(27r'v/— It)) is meromorphic on C, periodic with period 1, and quasi-periodic 
with period \og{p) /2tt^/ —\\ in other words, it is a theta function in the usual sense.) The canonical example of 
such a function is Op(z), a holomorphic p-theta function with multiplier —z^^; indeed, any holomorphic p-theta 
function can be written as a product of functions 9p{uz), and any meromorphic p-theta function as a ratio 
of such products. In the special case of multiplier 1, the function is called p-elliptic, for similar reasons. By 
standard convention, a p-theta function, if not explicitly allowed to be meromorphic, is holomorphic; however, 
p-elliptic functions are always allowed to be meromorphic (since a holomorphic p-elliptic function is constant). 
Let q be another complex number with \q\ < 1, such that p^ O q^ ~%. 

Definition 1. A p-theta q-dijference equation of multiplier /i(z) = az'' is an equation of the form 

viqz) = A{z)v{z), (2.1) 

where A{z) is a nonsingular meromorphic matrix (a square matrix, each coefficient of which is meromorphic on 
C*, and the determinant of which is not identically 0), called the shift matrix of the equation, such that 

A{pz) - ^l{z)A{z), (2.2) 

so in particular the coefficients of A are meromorphic p-theta functions of multiplier n{z). Similarly, a p-elliptic 
g-difference equation is a p-theta g-difference equation of multiplier 1. 

We will refer to the dimension of the matrix A as the order of the corresponding difference equation. We 
note the following fact about nonsingular meromorphic matrices. 

Proposition 2.1. Let M{z) he a nonsingular meromorphic matrix. Then M{z)^^ is also a nonsingular mero- 
morphic matrix, and if the coefficients of M{z) are meromorphic p-theta functions of multiplier ij,{z) , then those 
of M(z)^^ are meromorphic p-theta functions of multiplier ^{z)^^ . 

Proof. Indeed, the coefficients of the adjoint matrix det(M(z))M(z)~^ are minors of M(z), and thus, as polyno- 
mials in meromorphic functions, are meromorphic; this continues to hold after multiplying by the meromorphic 
function det(M(z))^^. For the second claim, if 

M{pz) = n{z)M{z), (2.3) 

then 

M{pz)-^ = ^jl{z)~^M{z)-^. (2.4) 

D 

Definition 2. Let v{qz) = A(z)v{z) be a p-theta g-difference equation. A meromorphic fundamental matrix 
for this equation is a nonsingular meromorphic matrix M{z) satisfying 

M{qz) = A{z)M{z). (2.5) 



It follows from a theoreiii of Praagman |15[ Thm. 3] that for any nonsingular meromorphic matrix A{z), 
there exists a nonsingular meromorphic matrix M{z) satisfying M{qz) — A{z)M{z) (this is the special case 
of the theorem in which the discontinuous group acting on CP is that generated by multiplication by q). In 
particular, any p-theta q-difference equation admits a meromorphic fundamental matrix. In the case of a first 
order equation, we can explicitly construct such a matrix. 

Proposition 2.2. Any first order p-theta q-dijference equation admits a meromorphic fundamental matrix. 

Proof. For any nonzero meromorphic p-theta function a(z), we need to construct a nonzero meromorphic 
function f{z) such that 

f{qz) = a{z)f[z). (2.6) 

Since a{z) can be factored into functions 9p{uz), it suffices to consider the case a(z) = 9p{uz), with meromorphic 
solution 

f{z)=Tj,,giuz); (2.7) 

this includes the case a{z) — hz^ by writing 

, _ ep{-hz)ep{-pzf-^ 
- ep{-bpz)ep{-zY-^- ^^-^^ 

D 

We note in particular that, since the elliptic Gamma function is symmetrical in p and q, the solution thus 
obtained for a first order p-theta q-difference equation also satisfies a g-theta p-difference equation. This is quite 
typical, and in fact we have the following result. 

Lemma 2.3. Let v{qz) — A{z)v{z) he a p-theta q-dijjerence equation of multiplier ^{z), and let M[z) be a 
meromorphic fundamental matrix for this equation. Then there exists a (unique) q-theta p-difference equation 
of multiplier fj,{z) for which M{zy is a fundamental matrix. 

Proof. An equation w{pz) = C{z)w{z) with fundamental matrix M{zY satisfies 

M{pzY = C{z)M{zY, (2.9) 

and thus, since M(z) is nonsingular, we can compute 

C{z) = M{pzYM{z)-K (2.10) 

This matrix is meromorphic, and satisfies 

C{qz) ^ M{pqzYM{qz)-^ = M{pz)A{qzY A{z)-^M{zy^ = Kz)C{z). (2.11) 



D 



By symmetry, we obtain the following result. 
Theorem 2.4. Let M[z) be a nonsingular meromorphic matrix. Then the following are equivalent: 



• (1) M{z) is a meromorphic fundamental matrix for some p-theta q-difference equation. 

• (2) M(z)* is a meromorphic fundamental matrix for some q-theta p-difference equation. 

• (1 ') M{z)^* is a meromorphic fundamental matrix for some p-theta q-difference equation. 

• (2') M{z)^^ is a meromorphic fundamental matrix for some q-theta p-difference equation, 

as are the corresponding statements with "some" replaced by "a unique". Furthermore, if the above conditions 
hold, the multipliers of the difference equations of (1) and (2) agree, and are inverse to those of (V) and (2'). 

Remark. In the elliptic case, the above observations were made by Etingof |8j, who also noted that the associated 
g-elliptic p-difference equation can be thought of as the monodromy of M. 

Given a p-theta q-difference equation, the corresponding meromorphic fundamental matrix is by no means 
unique, and thus we obtain a whole family of related (/-theta p-difference equations. There is, however, a natural 
equivalence relation on g-theta p-difference equations such that any p-theta g-difference equation gives rise to a 
well-defined equivalence class. First, we need to understand the extent to which the fundamental matrix fails 
to be unique. 

Lemma 2.5. Let M{z) and M'(z) be fundamental matrices for the same p-theta q-difference equation v{qz) = 
A{z)v{z). Then 

M'{z) = M{z)D{zf (2.12) 

for some nonsingular meromorphic matrix D{z) with q-elliptic coefficients. 

Proof. Certainly, there is a unique meromorphic matrix D{z) with M'{z) — M{z)D{zY , and comparing deter- 
minants shows it to be nonsingular. It thus remains to show that D{z) has g-elliptic coefficients, or equivalently 
that D{qz) = D{z). As in the proof of Lemma [^31 we can write 

A{z) = M{qz)M{zy^ = M'{qz)M'{zy^, (2.13) 

and thus 

D{qzYD{z)~* = M{qz)-'^M'{qz)M'{z)-'^M{z) = 1, (2.14) 

as required. D 

Theorem 2.6. Define an equivalence relation on q-theta p-difference equations by saying 

[v{pz) = C{z)v{z)\ ^ [v{pz) = C'{z)v{z)] (2.15) 

iff there exists a nonsingular q-elliptic matrix D{z) such that 

C'{z)D{z) = D{pz)C{z). (2.16) 

Then the set of q-theta p-difference equations associated to a given p-theta q-difference equation is an equivalence 
class. 



Proof. Let M{z) be a meromorphic fundamental matrix for the p-theta g-difference equation v{qz) — A{z)v{z), 
and associated g-theta p-difference equation w{pz) — C'{z)w{z). If Af'(z) is another meromorphic fundamental 
matrix for the g-difference equation, with associated p-difFerence equation w{pz) = C'{z)w{z), then 

Af (z) = M{z)D{z)\ (2.17) 

and thus 

C"{z) = M' {pzf M' {z)-^ = D{pz)M{pzfM{z)-^D{z)-'^ = D{pz)C{z)D{z)-\ (2.18) 

Conversely, if 

C'{z)D{z) = D{pz)C{z), (2.19) 

then M{z)D{zY is a fundamental matrix for the g-difFerence equation with associated p-difFerence equation 

w{pz) = c'{z)w{z). n 

Definition 3. The weak monodromy of a p-theta g-difFerence equation is the associated equivalence class of 
g-theta p-difFerence equations. Two p-theta g-difFerence equations are isomonodromic ii they have the same 
weak monodromy. 

Theorem 2.7. The p-theta q-difference equations v{qz) = A{z)v{z), v{qz) — A'{z)v{z) are isomonodromic iff 
there exists a nonsingular p- elliptic matrix B{z) such that 

A'{z)B{z) = B{qz)A{z). (2.20) 

Proof. Choose a g-theta p-difFerence equation w{qz) = C{z)w{z) representing the weak monodromy of the first 
equation. The equations are isomonodromic iff w{qz) = C{z)w{z) represents the weak monodromy of the second 
equation, iff the two equations have fundamental matrices satisfying M(qz) = M{z)C{z). But by Theorem 12.61 
(swapping p and q), this holds iff A'{z)B{z) — B{qz)A{z) For some p-elliptic matrix B{z). D 

Remark. Compare [5 , where the analogous result is proved for difference equations, relative to BirkhofF's [1] 
notion oF monodromy. 

Corollary 2.8. The map from isomonodromy classes of p-theta q-difference equations to their weak mon- 
odromies is well-defined, and inverse to the map from isomonodromy classes of q-theta p-difference equations 
to their weak monodromies. 

Note that the isomonodromy equivalence relation is also quite natural From the perspective oF the general 
theory oF difference equations (see, e.g. , [16| ) ; to be precise, two p-theta g-difFerence equations are isomonodromic 
iff they induce isomorphic difference modules. The latter fact induces a natural isomorphism between their 
difference Galois groups (assuming the latter are defined, effectively restricting this to the elliptic case). (Indeed, 
this preservation of Galois groups seems to be what is truly intended by the word "isomonodromy" , even in 
the differential setting, e.g., for non-Fuchsian equations, where the monodromy group conveys relatively little 
information.) 

A natural question, therefore, is the precise nature of the relation between the weak monodromy and the 
Galois group of a p-elliptic g-difference equation. For instance, both the isomonodromy class of the equation and 



its weak monodromy have well-defined Galois groups, and it can be shown (Etingof, personal communication) 
that in fact the groups are naturally isomorphic, with dual associated representations. (Thus, for instance, the 
fact that the difference equations we will be considering have triangular weak monodromy implies that they 
have solvable Galois group.) 

It is also natural to ask for a stronger notion of monodromy; for rational g-difference equations with suffi- 
ciently nice singularities, there is a well-defined notion of monodromy, an associated nonsingular g-elliptic matrix 
the nonsingular values of which generate a Zariski dense subgroup of the Galois group ([5]; see also Chapter 12 
of [16j). Krichever |13j defines an analogous matrix for generic difference equations with theta function coeffi- 
cients (although the relation to the Galois group is again unclear); although Krichever's genericity assumptions 
explicitly exclude the situation we consider above, his monodromy is again a difference equation with theta 
function coefficients. This suggests that the rational g-difference notion of monodromy should correspond at 
the elliptic level to a representative of our weak monodromy. In particular, this suggests the question of whether 
given a p-elliptic q-difference equation, there exists a representative of its weak monodromy the nonsingular 
values of which are Zariski dense in its Galois group. Another natural question is whether there is anything 
analogous to Krichever's local monodromy in our setting. 

It will be convenient in the sequel to introduce a slightly weaker equivalence relation. 

Definition 4. Two p-theta g-difference equations are theta-isomonodromic if there exists a nonsingular mero- 
morphic p-thcta matrix B(z) such that the shift matrices A{z), A'{z) of the equations satisfy A'{z)B{z) = 
B{qz)A{z). 

Theorem 2.9. Two p-theta q-difference equations are theta-isomonodromic iff their weak monodromies agree 
up to multiplication of the shift matrix by a factor of the form az^ . 

Note that v{qz) — A{z)v{z) and v{qz) — A{qz)v{z) are theta-isomonodromic with B{z) = A{z). 

3 Semiclassical biorthogonal elliptic functions 

In [21j . Spiridonov constructed a family of elliptic hypergeometric functions biorthogonal with respect to the 
density of the elliptic beta integral: 

{p;p){q;q) f no<r<6^,'7("'-^ ^) dz n -n r \ /qi\ 

o / — " r f.±2\ o r-r = 11 ^pA^rUs), (3.1) 

where the parameters satisfy the balancing condition 

Y[ Ur^pq, (3.2) 

0<r<6 

and the (possibly disconnected, but closed) contour is chosen to be symmetrical under z t— > 1/z, and to contain 
all points of the form p^q^Ur, i, j > 0, < r < 6, or more precisely, all poles of the integrand of that form. 

If we view this as the "classical" case, then this suggests, by analogy with [Ml [11] that we should study 
biorthogonal functions with respect to the more general density 

A^ ^(z;uo,...,U2m+5) == = — p , ^2\ ' (3-3) 



with new balancing condition 

n Ur^{pqr+\ (3.4) 

0<r<2m+6 

and the corresponding contour condition, integrated against the differential 



2ttJ^z 



(3.5) 



Note that if U2m+iU2m+b = M, then the corresponding factors of the density cancel, and thus we reduce to 
the order m — I density. Also, it will be convenient to multiply the integrands by theta functions, not elliptic 
functions, which has the effect of shifting the balancing condition. (The extent of the required shift can be 
determined via the observation that multiplying a parameter by q multiplies the integrand by a p-theta function; 
in any event, we will give the explicit balancing condition for each of the integrals appearing below.) 

One natural multivariate analogue of the elliptic beta integral is the elliptic Selberg integral [3 [TH] , the 
higher-order version of which we define as follows. 



iP,P)"{q;q)" f TT rp,,(iz±lz±l) ^ Ilo<r<2m+6U^4^) dz. 



.Am) ,^ ^, . ._ (P,p)"(g;g)" ( TT ^P,S4'4') n 



r,,,(t)-"2»n! Ja. ^J^^^ T.^zf^zf) ^^^\^ ^,,(41 2.^-lz. 

(3.6) 
where the parameters satisfy the conditions |t|, \p\, \q\ < 1, and 

i2"-2 -Q y^^(pq)-+i, (3.7) 

0<r<2m+6 

and the contour C is chosen so that C = C~^ , and such that the interior of C contains every contour of the form 
p^qHC, i,j > 0, and every point of the iormp^q^Ur, i,j>0,0<r< 2m + 6. (The latter set of points represents 
poles of the integrand; if (as often occurs below) some of these points are not poles, then the corresponding 
contour condition can of course be removed. Similarly, if the cross terms are holomorphic (e.g., if f = q, as is 
the case below), then C need not contain the contours p^qHC.) Note that if \uo\,. . . ,|u2m+5| < 1, then C can be 
chosen to be the unit circle. More generally, such a contour exists as long as p^qH^UrUs is never 1 for i, j, fc > 0, 
< r, s < 2m + 6, and the result is a mcromorphic function on the parameter domain. 
When m = 0, the elliptic Selberg integral can be explicitly evaluated: 

lln!t;p,q{uo,Ui,U2,U3,Ui,U5)^ ^ IJ,,g(f+^) Jl^ IJ,,g (f W^Ws) , (3.8) 

0<i<n 0<r<s<6 

while the order 1 elliptic Selberg integral satisfies a transformation law with respect to the Weyl group £'7; 
more precisely, the renormalized (holomorphic) function 

^^S;p,g(«o,...,W7):=//i;l;p,,(i'/'wo,...,i'/'u7) J] r+,.*(iu.w«) (3.9) 

0<r<s<8 

is invariant under the natural action of Ej on the torus of parameters |181 Cor. 9.11]. More importantly for our 
present purposes, when t ^ q, the renormalized order 1 elliptic Selberg integral satisfies an _E8-invariant family 
of nonlinear difference equations making it a tau function for the elliptic Painleve equation [17j Thm. 5.1] (for 



the relevant definition of tau functions, see [12] )• As an aside, it should be noted that TT] also showed that 
when t — q^/'^ or t — q^ , the integral satisfies slightly more complicated analogues of the tau function identities; 
as yet, neither a geometric nor an isomonodromy interpretation of those identities is known. 

Since we will be fixing p, q, and t = q in the sequel, we omit these parameters from the notation; we will 
also generally omit m, as it can be determined by counting the arguments. 

Consider the following instance of the elliptic Selberg integral: 

Fn{x;v) = x^'^v'^IIniuo, . . . ,U2m.+5,qx,pq/x,v,p/v), (3.10) 

satisfying, as usual, the balancing condition 

?'"-' n Ur^{pqr+\ (3.11) 

0<r<2m+6 

Since 






ipp{x,z):=x rp^q{qxz ,pq/xz ) = — p — -, — XiT ^ ^ 9p{xz ), (3-12) 



we see that the integrand of Fn{x;v) is holomorphic in x; indeed, it differs from the order m elliptic Selberg 
integrand by a factor 

In particular, the 2;-dependent conditions on the contour are irrelevant, as there are no a;-dependent poles. We 
thus find that Fn{x]v) is a i?Ci -symmetric theta function of degree n; that is, it is a holomorphic function of 
X satisfying 

F„{l/x;v)^Fr,ix;v) (3.14) 

F„ipx; v) = {px^y"F„ix; v) (3.15) 

(In general, BCn denotes the "hyperoctahedral" group of signed permutations, which will act by permutations 
and taking reciprocals.) This function satisfies a form of biorthogonality; to be precise, we have the following. 

Theorem 3.1. Let Gn{x) be any BC'i- symmetric theta function of degree n, and let C he any contour satisfying 
the constraints corresponding to the parameters uq, . . . , U2m+5, v,p/v with 

g2n-2 -Q Ur^{pqr+\ (3.16) 

0<r<2m+6 

Then for any x such that the contour C contains p^x and p*+^/x for all i > 0, 
{p-pf f Fn{z;v)Gn{z) ^,^.^, ^, , dz _^ ,^a^»+i..»+i ; 

(3.17) 



' It \l t N A(z;Mo,...,U2m+5)- 1==^ = Gn{x)x'^^ v'^+ IIn+l{U(i,...,U2m+5,X,p/x,V,p/v). 

c ipp(x,z)'ifjp{v,z) 27rv-lz 



In particular, i/iJ„_i(x) is a BCi- symmetric theta function of degree n~\, then 



[ Fn{z;v)Hn-i{z) dz 

/ —, ^ A(z;Mo,---,U2m+5)7; — P^^=0- (3-18) 



10 



Proof. If replace Fn{z;v) by its definition, the result is an n + l-dimensional contour integral over C"^^. 
Moreover, the integrand is very nearly symmetric between z and the remaining n integration variables. To be 
precise, if we write the original integration variable as zq, then the resulting integrand is a _BC„+i-symmetric 
factor multiplied by 

^"('°) (3.19) 



which is invariant under the subgroup BCi x BCn- If we average the integrand over BCn+i, this will not 

change the integral, as the contour is i3C„+i-invariant. We can thus replace the above factor by the average 

over cosets: 

1 v^ Gn{Zk) _ 1 Gn{x) 



the identity follows from the fact that if we multiply both sides by no<i<n ^p(^' -^«)i then both sides are 
-BCi-symmetric theta functions of degree n in x, and agree at the n + 1 distinct pairs of points z^ . 

The claim follows immediately. D 

Remark. At the level of orthogonal polynomials, such an n-dimensional integral representation is implicit in [24j : 
more precisely, Szego gives a representation of orthogonal polynomials as a determinant, but the Cauchy-Binet 
identity allows one to turn it into an n-dimensional integral involving the square of a Vandermonde determinant. 

Note that in the above calculation, the x-dependent constraint on the contour was only relevant to the 
eventual identification of the n + l-dimensional integral as an elliptic Selberg integral. We also observe that if 
V has the form Ur/ q, then the parameters Ur and p/v in the Selberg integrals multiply to pq and thus cancel. 
We thus find that Fn{z;Ur/q) satisfies biorthogonality with respect to a general order m instance of A(z). It 
will, however, be convenient to allow general v in the sequel. 

We thus see that the integral Fn{z; v) is in some sense an analogue of an orthogonal polynomial. Similarly, 
the n + l-dimensional integral of Theorem 13.11 is analogous to a Cauchy transform of Fn{z; v), as the integral 
of Fn{z;v) against a function with a moving pole. This suggests that these two integrals should form a row in 
the fundamental matrix of our difference equation. This leads to the question of how this row depends on v. 
Define 

F.;^{x,v) = v"+^x'^+'^^pp{v,x)IIn+liuo, . ■ .,U2m+5,X,p/x,V,p/v), (3-21) 

where the factor x"^^ is chosen to make the integrand invariant under x h^ l/x, the factor ii"+^ for symmetry, 
and the factor "tpp^v, x) to simplify the following identity. 

Lemma 3.2. The functions Fn{x;v) and Fj^(x;v) satisfy the identity 

Fn{x; v)F+{x, w) - Fn{x; w)F+{x, v) = //„(uo, . . . , U2m+z)F+{v, w) (3.22) 

Proof. Taking G„(x) = Fn{x;w) above gives 

F^{x;w)F:{x,v)^^-P^ f ^P(^;^) FUz;v)F^{z;w)A{z;uo,...,U2^+,) '^^ (3.23) 

2 Jc Wp(x,z)'ipp[v,z) 2tt^/-1z 

Thus the two terms on the left-hand side agree except in the first factors of the integrands; the difference of the 

two integrals can be simplified using the addition law, and gives a result independent of x; setting x = v gives 

the desired result. D 
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Similarly, we have the following. Let 

F~{x,v) := ipp{v,x)x^~'^v^~'^II„~iiuo, ■ . . ,U2m+5,qx,pq/x,qv,pq/v). (3.24) 

Lemma 3.3. For any BCi- symmetric theta function Gn of degree n, 

' A[z;uo,...,U2m+5)- — ■^=^ =Gn[x)Fn[v;x) -Gn[y)Fn[v;y). (3.25) 



2 Jc ^pix,z)i)p{y,z) ' '■'■' 27rV^2 

In particular, 

F~{x,v)F+{x,w) +Fn{v;x)Fn{x;w) = IIn{uQ, . . . ,U2m+z)Fn{v;w), (3.26) 

and if Hn-2{z) is any BGi-symmetric theta function of degree n — 2, 

dz 

F„ {z,v)Hn^2{z)A{z;UQ, . . . ,U2rn+5)- — p=^ = 0. (3.27) 

c 2n\' — lz 



(3.28) 



Remark. In particular, we see that 

ippiv,z) 
is essentially a biorthogonal function of degree n — 1. 

Theorem 3.4. The functions F^ , Fn and F^ satisfy the identities 

j_, ^ I Fn(x;u)\ ,, f Fn(x;v)\ , , , I Fn(x;w)\ 

P^i-^ ^) " - F^tiu, w) ; ' + f +(u, V) " U (3.29) 

\F+{x, u)J \F+{x, v)J \F+{x, w)J 

-^, N / Fl'{x, u)\ , ^ I F„(x;v)\ , , / Fn(x; w)\ , ^ 

F+{v, w)[ ; V M - F^,{u; w) "^ ' M + F„(u; v) "^ M - (3.30) 

\-Fn{u;x)) \F+{x,v)J \F+{x,w)J 

X / F~(x, u) \ , , / Fz^ix, v) \ , , F„ (x\ w)\ , , 

Fn{v; w)[ l\' \]~ F^iu; w) [ l\' \]+ F-{u, v) J) M = (3.31) 

\-Fn{u;x)J \-Fn{v;x)J \F+{x,w)J 

f F-(x,u)\ , ( F~(x,v)\ , { F-(x,w)\ 

Proof Each identity is the Phicker relation for the 2x3 matrix formed by concatenating the three column 
vectors that appear. In the first two cases, we have already computed the requisite minors; the remaining minor 
follows as a special case of the third identity, which can be derived by eliminating a common term from two 
instances of the second identity. D 

Remark. Note that the proof of these identities didn't require the balancing condition, or even that the biorthog- 
onality density was A. Furthermore, the only way in which the proof depended on properties of elliptic functions 
was in the fact that ipp satisfies a partial fraction decomposition result. If we generalize the results with this in 
mind, we find that these are precisely the generalized Fay identities of p^ll7j. 

We also note that the change of basis from F„ to F~ can be interpreted as relating degree n biorthogonal 
functions to degree n — 1 biorthogonal functions; i.e., the analogue of the three-term recurrence for orthogonal 
polynomials. 
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We thus sec that, as functions of x, the vectors 

( Fn{x;v) 



(3.33) 

and 

" ^ ' M 3.34 

\-Fn{v-x)J 

for aU w G C*, together span only a 2-diinensional space, and the change of basis matrix between any two such 

bases of this 2-dimensional space is computable in terms of F„, F,f . And, naturally, the choice of basis will 

have no effect on the resulting difference equation beyond conjugation by a matrix independent of x. Since it 

will be useful to allow an arbitrary such choice, we extend the notation by defining values for F,^ on hatted 

arguments (equivalently, defining F+ as a function on (C* W C*)^), as follows: 

F+{v,w):=F^{v;w), (3.35) 

F+{v,w):^-Fn{w;v), (3.36) 

F+{v,w):^F-{v,w), (3.37) 

F„{x;v):^F-{x,v)- (3.38) 

note that this extension of i^+ preserves its antisymmetry. Note that in this notation, the identities relating 
F„, F^ reduce to the single identity 

F+{w,x)F+{y,z)-F+{w,y)F+{x,z)+F+{w,z)F+{x,y)^0, (3.39) 

ior all w,x,y,ze C*l±)C*. 

To proceed further, we will need to understand how our functions behave under the monodromy action 
X t-^ px; it will also turn out to be useful to know how x i— > 1/x acts. Easiest of all is x i-^ p/x; in that case, the 
elliptic Selberg integral itself is manifestly invariant, so we need simply consider how the prefactors transform: 

F„{p/x; v) =. (xVp)"F„(x; v) (3.40) 

F+ip/x,v) - {p/x^rF+ix,v). (3.41) 

For X H^ 1/x, we similarly have 

F„{l/x;v)^F„{x;v). (3.42) 

However, for F^{l/x,v), while the integrand remains constant, the constraints on the contour change. Assume 

for the moment that v £ C* , and choose a i3Ci-symmetric theta function G'„(z) of degree n such that G„(a;) ^ 0, 

so that 

Fn{x,v)^ / — --— A z;no,---,M2m+5)- — 1==^- 3.43) 

2 Jc'4^p[x,z)ipp{v,z)Gn{x) 27rV^z 

Then x h^ 1/x leaves the integrand the same, but moves the contour through x and 1/x. Thus F^{l/x,v) — 

F^(x,v) can be computed by residue calculus; by symmetry, we find that it is twice the residue at z — 1/x: 

F+{l/x,v)-F+ix,v)^x-'9,{x^)Fn{x;v) [] Tp^giurX^'). (3.44) 

0<r<2m+6 

Putting this together, we obtain the following. 
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Lemma 3.5. The functions Fn and F^ have the monodromy action 

(F„(l/x;t;) F+il/x,v))^(^F„{x;v) F+{x,v))/l x-^e,ix^)ll,<r<2m+6^pA^rX^')\ (345) 

(f„(px;«) F+(px,i;)) = (f„(x;«) F+ix,v)) /(px^)-" iP^^r^-'Og{x^)Uo<r<2m+6^pA^rX^')\ 

\ (px2)" J 

(3.46) 
valid for all u G C* W C*. 

Proof. The only thing to check is that it extends to the other copy of C*, but this foUows immediately from 
the facts that the monodromy is independent of w G C* , and that for all u G C* W C* , the row vectors lie in the 
same 2-dimensional space. D 

This is not quite a g-theta p-difference equation as we would wish, but it is straightforward to turn it into 
a g-theta p-difference equation. Define a 2 x 2 meromorphic matrix M„(z; v, w) for w, w G C* W C*: 

, , / X ( Fn{z]v) z-i0p(z2)F+(z,i;)/A(z;uo,---,U2m+5)\ .„ .„. 

Mn\z\v,w) :— \ 19 • ('J-4') 

\Fn{z]w) Z^'^ep{z^)F+{z,w)/A{z;Ua,.. .,U2m+5)/ 

Theorem 3.6. The matrix M„(z; v, w) is a meromorphic fundamental matrix for ap-theta q-dijference equation 
with multiplier q^^". The isomonodromy class of the equation is independent of v and w, and invariant under 
all permutations of the parameters and all shifts 

(mo, • ■ • , M2m+5) ^ (g'^'UO, ■ • ■ , ?''"" + ■' M2m+5) (3.48) 

with kr (z li such that 

J2 kr = 0; (3.49) 

0<r<2m+6 

and invariant under simultaneous negation of all parameters. The theta-isomonodromy class is further invariant 
under all shifts 

(wo, . . . , U2m+5,z, n) H^ {q^°ua, . . . , q'"''"+''U2m+5, Q^z, n + v) (3.50) 

with I G 2^7 /c,- G Z + Z, V (zIj such that 

2v+ ^ kr^O. (3.51) 

0<r<2m+6 

In addition, the associated shift matrix A{z) satisfies the symmetry A{l/qz)A{z) ~ 1. 

Proof. Most of the claims follow immediately from the fact that M(z; w, w) satisfies the g-theta p-difference 
equation 

T,. I . ^, , .({pz^y" {pz^)"-^A{z;uo,...,U2m+5)/Mpz;ua,...,U2m+5)\ .„ ^„x 

Mn{pz;v,w)=M„{z;v,w)\ 2v.-2a. ^/A/ ^ ' ^^'^^^ 

V -{pz^y^ ^A(z;uo,. .. ,U2„i+5)/A[Pz;ua,.. .,U2m+5)J 
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and has determinant 

det(M„(z;u,u;)) = //„(mo, • •• , U2m+5)-F+(u, w)-r- -, (3.53) 

A(z;Uo,.. ■,U2m+b) 

so is nonsingular. The only additional thing to check for the isomonodromy claims is that 

(pZ^)^""^A(z; Mo, • • • , U2m+5) / ^{PZ; Mo, • • • , U2ni+b) (3.54) 

is invariant (up to scale, for theta-isomonodromy) under all of the stated transformations. The symmetry of 
A{z) follows immediately from the symmetry 

M(l/z; V, w) = M{z; v,w)i j . (3.55) 

D 

Remark 1 . One can avoid the appearance of theta-isomonodromy above at the cost of introducing some "ap- 
parent" singularities, and an additional parameter controlling the location of those singularities. Indeed, if one 
defines 

M^{z; x; v, w) := ^'1 — t7tM"„(z; v, w), (3.56) 

then M^ satisfies the same transformation law as M with respect to z i-^ 1/z, while under z h^ pz, one has 

M;(pz; x; V, m.) = M;(z; x; v,u.) (' ^^"Tl^f ^'' "°' ' ' ' ' "^"+^^/^^^"' ""'••" """+^^ ] . (3.57) 

\0 -{pz^y"^ ^A{z:Uo,...,U2m+5)/MpZ]Uo,...,U2m+5)/ 

Thus the associated shift matrix 

^-^'' ""' "' ") -^ ep{q^xz,x/qz)^-^'' "' ^) (^-^^^ 

is elliptic, with the same symmetry as A, and every shift 

(mo, • ■ • , M2m+5, z, n, a;) 1-^ ('7''°Mo, . . . , g''^'"+''M2m+5, 9'^, n + u^q'^ x) (3.59) 

with / G ^Z, Z', fcr e / + Z, z/ G Z such that 

2;y+ ^ fc^==0 (3.60) 

0<r<2m+6 

gives rise to a true isomonodromy transformation of this elliptic difference equation, with associated operator 
r,// N rp,,(/+'xz,/-'x/z) Tp,q{q''xz,q''xlz) 

B (Z\X\V,W) — -^—^ ; , , , ,, ; -77—; — -— Blz'.v ,w) . (3.61) 

^ ' ' ' ^ Vp^q{xz,x/z) i;,, (g^+^+W'xz, ?"+''+' -'x/z) ^ ' ' ^ ^ ' 

In particular, the isomonodromy transformations differ from the corresponding theta-isomonodromy transfor- 
mations by a meromorphic theta function factor depending only on v, I', I. There is also an isomonodromy 
transformation between ^^(z;x;w,w) and j4^(z;x'; w, w), for arbitrary x', but the corresponding B matrix is 
(generically) multiplication by an elliptic function of degree 2n. It follows that only those parameter shifts 
satisiying the integrality condition above can extend to arbitrary solutions of the elliptic Painleve equation (for 
which one effectively has noninteger n). 
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Remark 2. The lattice of possible kr vectors is the lattice £'2m+6 obtained by adjoining the vector 

(1/2, 1/2, . . . , 1/2, -1/2, -1/2, . . . , -1/2) (3.62) 

of sum to the root lattice -D2m+6- In particular, when m — 1, this lattice is precisely the root lattice Eg- 

4 The difference equation 

Naturally, simply knowing the existence of a difference equation with associated isomonodromy transformations 
is of strictly limited usefulness, so we would like to be more explicit about the equation, and at the very least 
generators of the group of monodromy-preserving transformations. 

The first thing we will need to understand about the shift matrix is the locations of its singularities; i.e., the 
points where the coefficients have poles or the determinant has a zero. This in turn depends on determining 
the divisor of F+. Define 

{x;p,q):= 1[{1 -p'q'x), (4.1) 

0<i,j 

with the usual multiple argument conventions. 
Lemma 4.1. The function 

Yl iurx,pur/x;p,q)\ F+{x,v) (4.2) 

y0<r<2m+6 / 

is holomorphic for x G C*. If v Cz C* , the function vanishes at x — v, x = p/v. 
Proof. As before, assuming w G C*, we have 

2 Jc i'p\XiZ)ipp[v,z)Gn{x) 2ii^/-lz 

(jp]pf f XV^pp{v,x)Fn{z;v)Gn{z) / \ d,^ fA A\ 

I '^ -A(z;'Uo,---,W2m+5,a;,p/a;,w,p/w)- — -j==^. (4.4) 



2 Jc Gn[X) 2'K\/-\Z 

where G„(z) is any _BCi-symmetricp-theta function of degree n not vanishing at x. Since Fn{z; v) is holomorphic 
in z, we may apply Lemma 10.4 of |18j (note that condition 3 of that lemma reduces in our case to the balancing 
condition) to conclude that 

JJ {urX,Urp/x;p,q) {xv,px/v,pv/x,p'^/xv;p,q){xv)~'^'ipp{v,xy^Gn{x)F+{x,v) (4.5) 

,0<r<2m+6 / 

is holomorphic in x. (The conclusion concerning the x-independent poles is not useful to us, as Fn{z] v) certainly 
has singularities that depend on the remaining parameters.) This nearly gives us the desired result, except for 
the factor G„(a;), which disappears by the fact that F^{x, v) is independent of G„, and the additional factor 

{qxv , pqx / V , pqv / X, p'^ q/ xv; p, q). (4.6) 

This latter factor can be eliminated, and the result extended to w G C* W C*, by expressing F^{z, v) as a linear 
combination of Fj^(z,w) and ^,^(2,1/;'), which for generic w and w' G C* are holomorphic at the offending 
points. n 
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Theorem 4.2. The matrix 

An{z;v,w) := {q-^z-^ Y[ Op{urz) \ An{z;v,w) (4.7) 



0<r<2m+6 



q'^z-^ n Op{urz)\ M„{qz;v,w)M„iz;v,w)-^ (4.8) 

0<r<2m+6 / 

is holomorphic, with determinant 

Aci{An{z;v,w))= ]^ 0p{urZ,Ur/qz), (4.9) 

0<r<2m+6 

satisfies the p-theta transformation law 

A^{pz; V, w) = {pqz^)-"'-^An{z; v, w), (4.10) 

and has the symmetry 

A„{l/qz;v,w)=l \An{z;v,wYi^ 1. (4.11) 

Proof. The formula for det { An {z;v,w)) follows immediately from the formula for dct(Mn{z; v , w)) . Similarly, 
the fact that An{l/qz;v,w)An{z;v,w) — 1 becomes 

An{l/qz;v,w)An{z;v,w) = det A„{z;v,w); (4-12) 

the symmetry of A„(z; v, w) follows from the usual formula for the inverse of a 2 x 2 matrix: 

C-=det(C)-f^;^ -':A (4.13) 

y — C21 Oil / 

Another use of this formula allows us to explicitly write down the inverse of Mn{z;v,w). This, in turn, 
allows us to express the entries of A{z) as polynomials in Fn and F,^ with coefhcients that are (holomorphic) 
p-theta functions in z: 

^^( )^/K(<Z.;.) FnHqz,v)\ Uz) oUF^tM -FnHz,v)\ 

\Fn{qz;w) F+{qz,w)l I b{z) \-Fn{z;w) Fn{z;v) I 



F,i{qz;v) b{z)F+{qz,v)\ fa{z)F+{z,w) -a(z)i^+(z, u) 
Fn{qz;w) b{z)F+{qz,w) I \ ~Fn{z;w) Fn{z;v) . 



(4.15) 



where 



g ^ llo<r<2m+6"p("T-^^) 
//„(mO, • • • , U2m+h)Fn{v, w) 



«(^) = TTT- ^. ^T.+ ( ^ (4.16) 



IIn[Uo, ■ ■■ , U2ni+5)Pn{v, W) 
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In particular, the only possible poles of An come from poles of a{z) F^^ (z , v) , a(z)_F'^(z, w), b{z)F^{qz,v) 
and b(z)Fjl{qz,w), or, equivalently, poles of 

a{z) n {urz,pur/z;p,q)-'^ n , (P/"-";P) (4.18) 

0<r<2m+6 0<r<2)n+6 ^ ^ ' -^ ' '-^'^^ 

and 

h{z) n K'z^,P^i./9^;p,9)-^- n , (g"/"-;P) (4.19) 

0<.<2™+6 0<.<2™+6 i^rqZ,pUrlz-p, q) 

It follows that 

W^ {urqz,pUr/z;p,q)An{z;v,w) (4.20) 

0<r<2m+6 

is holomorphic in z. But the entries of An{z;v,w) are meromorphic p-theta functions, and thus their divisors 
are periodic in p. Since the remaining set of potential poles contains no p-periodic subset, there are in fact no 
surviving poles. D 

We can also compute the value of An{z; v, w) at a number of points. 

Theorem 4.3. The matrix An{z;v,w) has the special values 

7, , ^ <lu7^Ila<r<2m+6^pi^rUs/q) /^f„(u,;i;)^ (F+{us/q,w) -F+{us/q,v)) 
An{us/q;v,w) ^ —- = — p- rl , ,p A (4.21) 

I ,-,1 ^ g^s''no<r<2r»+6^pK"^/g) /^ -F+ (u./g, ^;) ^ (-F„(m,;w) F„(m,;w)) 

//„(uo, • • . ,U2™+5)iV(w, w) \F+{us/q,w)J 

for < s < 2m + 6. In addition, we have the four values 

Aniq-'/^;v,w)= n 9piurq-'/^)r °) (4.23) 

0<r<2m+6 V / 

M-q-'/';v,w)= n 9,{-Urq-'^^) r °) (4.24) 

0<r<2m+6 V / 

in((p/<z)'/';«,ti') = ?->-' n (^piMp/qy^^)[ J (4.25) 

M-{p/qy/';v,w)^q-y-^ n ^p(-"r.(p/g)^/') r "j (4.26) 

0<r<2m+6 V^ / 

at the ramification points (fixed points of z i-^ ^/qz modulo (p)). 

Proof. We first observe that aX z — Ug/q, b{z)F^{qz,v) and b{z)F^{qz,w) vanish, and thus the formula for 
A{us/q) simplifies as stated. The second set of special values follows similarly from the vanishing of a(z)F+(z, v) 
and a{z)Fjl{z, w) at z — p/us, together with the p-theta law of A. 

When z = ztq^^''^, so that qz ~ 1/z, we find a{z) — b{z), Fn{qz,v) — Fn{z,v), and F^{qz,v) — F^{z,v); 
the last difference vanishes due to the factor Oq{z^) = in the relevant residue. The expression for A{±q^^^'^) 
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thus simplifies immediately. Similarly, at z = ±^/p/q, we have qz — p/ z, and again the entries immediately 
simplify. D 

Note that the symmetry of A and the elementary values at the ramification points imply that the matrix is 
already determined by its values at Us/q for any m-\-2 values of s (assuming u^ are generic); the above special 
values are thus highly overdetermined. It is also worth noting that if v and w are of the form Ur/q or u^ (with 
different values of r), then all but two pairs of special values can be expressed entirely in terms of order m 
elliptic Selberg integrals with shifted parameters. This specialization also has the effect of causing the kernel 
of An(ur/q) and image of An{l/ur) (or vice versa, as appropriate) to be coordinate vectors. 

We also note that when n — Q, Fo{z;w) = F^^{z,w) — 0, and thus if v is not "hatted", then Ao{z;v,w) is 
well-defined and triangular; note in particular that Fq{v,w) — 1. In particular, it follows that for any n > 0, 
An{z] V, w) is isomonodromic to a triangular shift operator with at most the same number of singularities. 

It remains to consider the isomonodromy transformations. Changing v and w is straightforward, as we have 
already observed; the precise isomonodromy transformation follows from 

Mn{z;v ,w)^-— r ,, , , ,, , ] M^{z;v,w). (4.27) 

With this in mind, we can feel free to make choices for v and w if this will simplify the expressions for the 
remaining isomonodromy transformations. 

We first consider the case of integer shifts, or in other words shifts under the lattice D2m+6- It thus suffices 
to consider the two cases {uq, mi, n) i-^ (qwo, ui/q, n), {quo, qui, n — 1). 

Lemma 4.4. We have the (isomonodromy) transformations 



— /I \ _ 

M„(z; Uq, Ui/q; QUq, Ui/q, U2, . . . , M2m+5) == ("O?/-"!)" I 0p(«iz±V9) I ^'^niz; Ui/q, Uo; Uo,..., U2m+5) (4.28) 





Mn^i{z;uo,Ui;quo,qui,U2,...,U2m+5) ^ iuQUiY' i ''p("i^ ) ^^ J Af„(z; uj, uj^; Uo, . . . , M2m+5) 

(4.29) 

Remark. In the version with apparent singularities controlled by x, the transformations (j4.27p and (j4.28p remain 
unchanged, while (I4.29P should be multiplied by 9p{q"~-^xz^^). 

It remains to consider the case of a half-integer shift, say 

{uq, . . . , Um+2,Um+3, . . . , W2m+5, Z, V, w) 

^ (ql/^uo, . . . , g'/'um+2, q-^^^Ura+3, ■■■, q-^'''u2rn+5, q^'^ z , q^'^v', q^'^w'). (4.30) 

In this case, just as with A itself, we cannot give a closed-form expression, so must be content with a description 
of the singularities and a sufficiently large set of special values. Let Bn{z] v, w; v' , w') denote this isomonodromy 
transformation; that is, 

Bn{z; V, w; v' , w') = Mn{q^/'^z; q^/'^v', q^/'^w'^u'^, ..., U2„+5)M„(z; v, w;uq,..., U2m+5)"\ (4.31) 
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where 

{q^^^Ur 0<r<TO + 3 
(4.32) 
q^^/'^Ur m + 3 < r < 2ra + 6. 

Note that 

Bn{l/qz; V, w] v\ w'j^^Bniz; v, w; v' , w') — An{z; v, w). (4.33) 

For convenience, we define 

Gtiv, w) :- F+{q'/^v, q^'^w; u'„, . . . , u',^^,), (4.34) 

and similarly for G'„(u; w); this will free us to again omit most of the parameters from Fn and F^ . 
Theorem 4.5. The matrix 

B„{z; V, w; v', w') := I (q'^^z)-' JJ Opiurz) B„{z; v, w; v' , w') (4.35) 

y 0<r<m+3 J 

is holomorphic in z with p-theta law 

B^ipz; V, w; v' , w') = I q-y-\-ir+^ J] V^ ] z-"-3^„(z; v, w; v' , w'), (4.36) 

has determinant 

-q ' -pr-, TTTT r 11 ^pi^rZ) [[ 0p{Ur/qz), (4.37) 

and has the special values 

f,,-,, , ,s -^i^Ilm+3<r<2m+e^piUrUs/q) ( G+{Us / q,v')\ {-Fn{u,-w) F„(u,;«)) 

//„(U0, • • • , U2rn+b)Fn(v, w) \G+{Us/q, w')J 

for < s < m + 3, and 

P, , , ,, (g-^/^^,)-l^o<.<n^+3^pK^^/g) fGn{us/q;v')\ (F+{uJq,w) -F+{u,/q,v)) 
Bn(us/q;v,w-v ,w)^ — = — — ^ / ,xP ^ 

IIn(Uo,.. .,U2,n+5)F,T{v,W) \Gn{Us / q;w') ) 

(4.39) 
for ?Ti + 3 < s < 2m + 6. 

Proof. As in the computation for A„, we can use the known determinant of A'/„ to write the entries of _B„ in 
terms of F„ , F^ , G„ , and 6*+ . We find 

~ , (G^{z-v') d{z)Gt{z,v')\ (c{z)F+{z,w) -c{z)F+{z,v)\ 

Bn{z;v,w,v ,w = , 4.40) 

\Gniz;w') d{z)G+{z,w') \ -F^(z;w) F^{z;v) /' 



where 

iq^^^z)-^llo<r<m+3^piUrZ) 
IIn{uQ,. ..,U2ra+5)Fn{v,w) 



c{z) = Z , ^^"<'-<"'+;j^;^^-"; ^ (4.41) 



d[z) = ^ nm+3<r<2m+6 SpJUr/qz) ^ ^^^^^^ 

IIn{uQ, . . .,U2m+5)Fn{v,w) 
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Since 

C{Z)F+{Z V) Uo<r<,n+3iP/^rZ;p) ^^ ^^^ 

d{z)GUz.v') ^ = -. n^+3<.<.^+a(g^/^.;P) (4 44) 

\A^<r<,n+-,ilUrZ,pUr/z-p, q) \Am+:i<r<2m+QV^^rZ,pUr / Z'p, q) 

and the coefficients of Bn are p-tlieta functions, we conclude as before that _B„ is holoniorphic. 

The special values again follow by choosing z so that c{z) or d{z) vanishes, and using the p-theta law as 
appropriate. D 

Remark 1. The relation between A and B becomes, via the usual expression for the inverse, the expression 

I B„(l/qz;w,w;w',w')* j i?„(z; u, u;; u', w') = Ci„(z; w, w), (4.45) 

where 

p ^ -q-^'^IIn{u'o,...,u'^^+^)GtW,w') (4 4g) 

IIn{uo,.. .,U2m+5)Fn{v,w) 

Remark 2. In the form with apparent singularities, i?„ gets multiplied by 6p{xz)/0p{q"xz), making it elliptic, 
as expected. Moreover, since for each generator of our lattice of shifts we have exhibited an isomonodromy 
transformation with coefhcients of degree independent of n, the same holds for an arbitrary shift. 

Once again, the special values, together with the determinant and the p-theta law, are more than sufficient 
to determine -B„; indeed, for generic parameters, any to + 3 of the special values suffice. For both v4„ and i3„, 
this gives rise to a number of relations between the coefficients. For instance, using _B„, we can express 

G„ {Urn+s/q; V)F^ {Urn+s/q, w) (4.47) 

as an explicit linear combination of the terms 

G+ius/q,v')Fn{us;w) (4.48) 

for 0<s<7TT, + 3. If we choose v' and w suitably, we can arrange both for some terms of the resulting identity 
to drop out, and for the remaining integrals to be order m elliptic Selberg integrals. For instance, if v' = uo/q, 
w = ui , then 

GtAuo/q,v') = Fn{ui;w) = 0, (4.49) 

and thus the s = 0, 1 terms disappear, and we are left with an expression for 

Gn{um+3/q;uo/q)Fn{ui;Ujn+3/q) (4.50) 

as a linear combination of 

G+{u,/q,uo/q)F-{u,,ui) (4.51) 

for 2 < s < "^ + 3. When m — I, these identities and the corresponding identities arising from the 12 entry 
of An give direct proofs (i.e., without using the E-/ symmetry) of new special cases of Theorem 5.1 of |17] 
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(which states that the ehiptic Selberg integral satisfies bihnear relations making it a tau function for the elliptic 
Painleve equation). 

The special case B{z;ue/q,U'//q;uo,ui) (with tti = 1) is particularly nice. In that case, in contrast to the 
situation with A, all of the integrals that appear in the expressions for the singular values are order 1 elliptic 
Selberg integrals, and may thus be expressed via (|3.9p as tau functions for elliptic Painleve. It should follow 
from a suitable Zariski density argument that the various formulas resulting from consistency of this expression 
and of the action of isomonodromy transformations continue to hold for arbitrary tau functions, with any 
appearance of 5" replaced by Q G C* such that 

Q'q-' n ^r^{pqr+'. (4.52) 

0<r<2m+6 

(Sketch: As n varies over (large) positive integers, the balancing condition describes a dense countable family of 
hypersurfaces in parameter space; that the contour integral is dense among all solutions on such hypersurfaces 
follows from the fact that its difference Galois group is generically equal to GL2.) However, since the arguments 
of [3] give a much more conceptual proof of this fact (and, conversely, that any function satisfying all consistency 
conditions is a tau function), there seems little point to fleshing out the details of the Zariski density argument. 
Also of interest in the case m = 1 is the relation to the action of the Weyl group E-j. This turns out to 
be easiest to describe in terms of the B matrices, although what meaning this has from a difference equation 
perspective is as yet unclear. 



Theorem 4.6. In the case m = 1, let 



whe 



q^^/'^Ur/x < r < 4 
q^/'^UrX 4 < r < 8, 

1/4 / _2-« \ 1/2 /„, „, „, „, X 1/2 



f uQUiU2U'i \ _ f pq^ " \ _ / 



U0U1U2U3 



\u4U5UeU7 J \U4U5UGUr J \ pq^ " 

Then there exist matrices C and D independent of z such that 



(4.53) 



(4.54) 



Bn{z\ V, w; v', w'; uq,. . . , uj) = CBn{q^^'^xz; v" , w";v"', w"';uq, . . . , u'^)D. (4.55) 

Proof. Indeed, it suffices to verify this in the special case 

V — Ue/q w — u-i/q v' ~ uo/q w' = Ui/q (4.56) 

v"^u'e/q w"^uyq v"'^uyq w'" ^ u[/q, (4.57) 

for z — Ur/q, 4 < r < 7, in which case it follows readily by several applications of the transformation law |181 
Thm. 9.7] 

//„(uo,...,M7)= n ( n ^>Ai"''urus) n rp,,(9"--'w,M,) //„(u^,...,u^). (4.58) 

l<i<" \0<r<s<4 4<r<s<8 J 

□ 
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Finally, we observe that the case m = is precisely the elliptic hypergeometric equation [23j for a (ter- 
minating) elliptic hypergeometric series, since in that case the biorthogonal functions are Spiridonov's elliptic 
hypergeometric biorthogonal functions. 
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